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Abstract 

The traditional Power Grid has been designed in a hierarchical fash- 
ion, with energy pushed from the large scale production factories towards 
the end users. With the increasing availability of micro and medium scale 
generating facilities, the situation is changing. Many end users can now 
produce energy and share it over the Power Grid. Of course, end users 
need incentives to do so and want to act in an open decentralized energy 
market. In the present work, we offer a novel analysis of the Medium 
and Low Voltage Power Grids of the North Netherlands using statistical 
tools from the Complex Network Analysis field. We use a weighted model 
based on actual Grid data and propose a set of statistical measures to 
evaluate the adequacy of the current infrastructure for a decentralized en- 
ergy market. Further, we use the insight gained by the analysis to propose 
parameters that tie the statistical topological measures to economic fac- 
tors that influence the attractiveness of participating in such decentralized 
energy market, thus identifying the important topological parameters to 
work on to facilitate such open decentralized markets. 



Keywords: 



Power Grid, Decentralized energy trading. Complex Network Analysis 



1 Introduction 

The Power Grid is one of the engineering masterpieces of the XIX-XX cen- 
tury, being one of the most important infrastructures that contributes to the 
economic growth and welfare of any country. It has been designed as a hierar- 
chical system with large generating facilities on top, and a pervasive network of 
cables to distribute the energy to the geographically dislocated end users. Tra- 
ditionally, it has been created to be managed by a monopolist or an oligarchy 
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of actors. Typically, energy availability is given for granted, though its impor- 
tance becomes well too apparent both at the household and country level when 
prolonged blackouts strike and electricity flow is interrupted [TJ [5] . 

Something is though changing in the way energy is produced and distributed 
due to both technological advancements and the introduction of new policies. A 
clear trend of market unbundling is emerging (e.g., [21131) entailing the addition 
of many new players to the energy sector with the possibility to produce, sell and 
distribute energy. From the technological perspective, new energy generation 
facilities (mainly based on renewable sources) are becoming widely accessible. 
These are convenient and available at both the industrial and the residential 
scale ini [6]. The term Smart Grid, which does not yet have a unique agreed 
definition, is sometimes used to define the new scenario of a Grid with a high 
degree of delocalization in the production and trading of energy. The new actors, 
who are both producers and consumers of energy, referred to as prosumers, are 
becoming more numerous and will most likely demand a market with total 
freedom to energy trading [7] . In this coming scenario the main role of the High 
Voltage Grid may change, while the Distribution Grid (i.e., Medium Voltage 
and Low Voltage end of the Power Grid) gain more and more importance, while 
requiring a major update. In fact, the energy interactions between prosumers 
will increase and most likely occur at a rather local level, therefore involving 
the Low and Medium Voltage Grids. 

Given such emerging scenario, we propose to look at the lower layers of 
the Power Grid in a statistical manner, considering global metrics coming from 
the field of Complex Network Analysis. Few such studies exist in the literature 
using unweighted models of the High Voltage Grids. These have been performed 
especially to establish the resilience to failures of critical national infrastructures. 
It will not surprise that often such studies have appeared immediately after a 
major national blackout. More recently, these techniques have also been used for 
generating tests for the Smart Grid in [5] . Here we propose a novel study of the 
properties of the Medium and Low Voltage networks using the Northern part of 
the Netherlands as data source. Our investigation goes beyond the study of the 
existing as it also proposes a way of evaluating the infrastructure in terms of its 
ability to support delocalized energy trading. We argue that global topological 
statistical measures do influence the eagerness of prosumers to trade energy, 
which in turn, entails a structural modification of the Power Grid. 

The paper is organized as follows. Motivation for the present study and an 
initial overview of the state of the art are presented in Section [2j In Section [3j 
we introduce the data set used for the study and provide essential definitions. In 
Section|3]we provide an initial unweighted Complex Network Analysis. Section[5] 
is dedicated to a weighted Complex Network Analysis which provides the best 
insight on the available sample and a set of novel measures in the statistical 
study of the Power Grid. A comparison with an unweighted Complex Network 
Analysis is provided Section |6l The proposal on how to tie topological metrics 
to economic factors is presented in Section [71 together with an example. A 
detailed account of related work on Complex Network Analysis for the (High 
Voltage) Power Grid is provided in Section [5| Final discussion and conclusions 
are offered in Section [Hi 
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2 A Statistical View on the Power Grid 



Many scientific knowledge areas contribute to the design and analysis of Power 
systems: Physics (eleetromagnetism, classical mechanics), Electrical engineer- 
ing (AC circuits and phasors, 3-phase networks, electrical systems control the- 
ory) and Mathematics (linear algebra, differential equations). Such traditional 
studies tend to have a "local" view of the Grid, e.g., defining how to design 
a transformer and predicting its functioning. Typically studies tend to focus 
on the physical and electrical properties (e.g., O), or the characteristics of the 
Power Grid as a complex dynamical system [lOj . or again, the control theory 
aspects [TT]. 

The move from a "local" to a "global" view of the Power Grid as a complex 
system is possible by resorting to statistical graph theory, better known as 
Complex Network Analysis ( CNA ), that is, a statistical analysis of the dynamics 
of large graphs with the goal of identifying their characterizing properties, such 
as the average path length between any two nodes. Studies of this sort have 
been recently performed for example on the American [121 US] : the Italian [M] , 
and the Scandinavian |15j Grids; the European Grid as a whole is analysed 
in [161 117] . Interestingly, these studies, though using the same mathematical 
machinery, reach different conclusions on the characteristics of the identified 
graph (e.g., node degree distribution, network category). An explanation for this 
is that different geographical infrastructures may indeed have different layouts, 
thus yielding different topological properties. 

Motivations for such CNAs are common to all previous researches and are 
mainly two: identifying the right complex network model for the Power Grid; 
and provide a resilience analysis for predicting blackouts and critical components 
of the infrastructure. Folkloristic is that Power Grid CNA studies are 'popular' 
after a major blackout occurs, such as the North American black-out of 200^13 
or the Itahan one of 200^ (e.g., [H [H [H [T2]). The fragility of the Power 
Grid has been the major reason of concern that has determined the focus of 
such CNA studies on the High Voltage. In fact. High Voltage failures impact 
a large part of the network, thus resulting in electricity service disruptions of 
large portions of a whole country. From a graph perspective, the High Voltage 
samples for an entire country are usually quite small in terms of order and 
size. For instance, the Grids connecting the 15 European countries analysed by 
Rosas-Casals et al. in [T7] are composed overall of about 2700 nodes and more 
than half of the analysed samples are below 100 nodes; moreover, the graphs 
studied in [191 \7U\ are well below 200 both nodes and edges, while the sample 
analysed in [M] has less than 400 nodes and a little more than 500 edges. 

The motivation for the current work is quite different from previous ap- 
proaches. We consider the Power Grid as an infrastructure for decentralized en- 
ergy exchange. To this end, we are interested in the properties of the Medium 
and Low Voltages. Furthermore, we consider the Power Grid not only in its 
basic topology to assess resilience and connectivity, but also in the physical 
properties of the lines to assess the capacity of the infrastructure in supporting 
distribution. To the best of our knowledge, the study has novel motivations, in 
line with the current trends of the Smart Grid, but also studies new networks 
with novel weighted models. 

^http : //news .bbc . co .uk/2/hi/emier lcas/3152451 ■ stmj 
^ihttp : //news . bbc . co . uk/2/hi/3146136 . stmj 
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3 Northern Netherlands Medium and Low Volt- 
age Complex Network Analysis 



We focus on the Medium and Low Voltage Power Grid networks of (Northern) 
Netherlands. This choice is dictated by the fact that the Netherlands has a 
modern infrastructure. The Dutch High Voltage Power Grid is owned and man- 
aged by one player, Tennet, while the lower layers are partitioned geographically 
among fourteen companicfo that have their own distribution network across the 
country. The partition of the territory among energy distribution companies is 
shown in Figure [TJ 



Figure 1: Distribution companies over the Netherlands. Each color corresponds 
to one company. (Source: www.energieleveranciers.nl) 

The Grid information used in this study is provided courtesy of Enexis B.V., 
the distribution operator of Northern Netherlands. The provided data includes 
information about the transformers in the Grid together with the distribution 
substations. The data set also provides information about the distribution lines 
used to connect substations containing the length of cable and other interesting 
physical properties (e.g., resistance, capacity, voltage). For the sake of precision, 
we define the notion of a (Weighted) Power Grid graph. 

• All the substations and transformers arc considered equal and are repre- 
sented as nodes of the graph. 

• The cables connecting the substations are considered equal despite the 
differences in voltages and current carried and their physical properties, 
and thus modelled as unweighted edges in the graph. 

■^Zonc identifier and distribution provider: 1) RENDO Netwcrken, 2) Cogas Infra cn Be- 
heer, 3) Liander (former Continuon Netbeheer), 6) Stedin (former Eneco), 7) Westland Infra, 
8) ONS Netbeheer (now Stedin), 9) DELTA Netwerkbedrijf, 12) NRE Netwerk, 13) Enexis 
(former Essent Netwerk), 14) InfraMosane (now Enexis). 
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Present study 


Random Graph 


ID 


Order 


Size 


Avg. 
d 


APL 


GPL 


7 


APL 


GPL 


7 


1 


17 


18 


2.118 


3.398 


3.313 


0.00000 


1.427 


1.688 


0.13726 


2 


15 


16 


2.133 


3.086 


3.000 


0.00000 


2.319 


2.358 


0.00000 


3 


24 


23 


2.087 


4.499 


4.228 


0.00000 


3.127 


3.091 


0.05508 


4 


30 


29 


1.933 


4.545 


4.449 


0.00000 


1.860 


2.242 


0.05778 


5 


188 


191 


2.032 


17.726 


17.878 


0.00000 


3.846 


4.345 


0.00532 


6 


10 


9 


1.800 


2.423 


2.223 


0.00000 


0.978 


1.167 


0.26667 


7 


63 


62 


1.968 


5.204 


5.404 


0.00000 


2.514 


2.904 


0.03175 


8 


28 


27 


1.929 


4.784 


5.000 


0.00000 


2.553 


2.945 


0.04762 


9 


133 


140 


2.105 


11.543 


11.366 


0.01112 


3.702 


4.172 


0.01482 


10 


124 


138 


2.226 


8.053 


7.070 


0.00869 


3.010 


3.540 


0.02914 


11 


31 


30 


1.935 


4.353 


4.357 


0.00000 


1.590 


1.969 


0.07475 



Table 1: Low Voltage samples from the northern Netherlands Power Grid com- 
pared with Random graphs of the same size. 

• For the data samples that present disconnected components, each compo- 
nent is treated distinct graph. 

• The edges are considered undirected. 

These assumptions are common in Power Grid analysis from a graph theoretic 
perspective, see for instance [IHl [HI [HI [IHl [HI [H] and lead to the following 
definition. 

Definition 1 (Power Grid graph). A Power Grid graph is a graph G{V,E) 
such that each element Vi &V is either a substation, transformer, or consuming 
unit of a physical Power Grid. There is an edge Cij = {vi,Vj) G E between two 
nodes if there is physical cable connecting directly the elements represented by 
Vi and Vj . 

The next step is to bring cable properties into the graph definition. 

• For each cable connecting elements in the Grid a weight is defined based 
on the multiplication of the following quantities: 

— The principal resistance characterizing the cable (whose value is given 
in Ohm/km). 

— The length of the cable (whose value is given in km). 

• A special kind of connection is defined in the Power Grid known as a 
'link'. These are connections, usually very short, with negligible resistance 
for which the specific value is not provided in the dataset. For edges 
representing these links a conventional weight of 10~^ is given. This docs 
not affect the overall validity of the weighted model since the number of 
links in a sample is extremely limited (about 1% of the overall connections 
are made of links). 

Definition 2 (Weighted Power Grid graph). A Weighted Power Grid graph is a 
Power Grid graph Gw(V,E) with an additional function / : — > R associating 
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a real number to an edge representing the resistance, expressed in Ohm, of the 
physical cable represented by the edge. 

The analysis we perform uses samples from the Low Voltage and Medium 
Voltage Grids. The Low Voltage samples sum up to a total of 663 nodes and a 
683 edges; while the Medium Voltage samples sum up to 4185 nodes and a 4574 
edges. The size of the data set, tough being a sample and not the whole network, 
is about the same size or larger than those used in other available studies on the 
(High Voltage) Power Grid [HJdllllllilllTllIlKin]. We begin our analysis 
by considering the unweighted model to derive basic topological properties and 
then proceed with a richer investigation by introducing graph weights. 
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Random Graph 


ID 


Order 


Size 


Avg. 

d 


APL 


GPL 


7 


APL 


GPL 


7 


1 


191 


207 


2.168 


9.288 


8.990 


0.00296 


4.616 


5.079 


0.00225 


2 


884 


1059 


2.396 


9.817 


9.527 


0.00494 


5.440 


6.010 


0.00170 


3 


444 


486 


2.189 


11.033 


10.858 


0.00537 


5.547 


6.163 


0.00333 


4 


472 


506 


2.144 


17.095 


17.174 


0.01360 


5.039 


5.700 


0.00106 


5 


238 


245 


2.059 


11.715 


11.580 


0.00000 


3.558 


4.234 


0.00595 


6 


263 


288 


2.190 


12.775 


12.311 


0.01118 


5.046 


5.368 


0.01080 


7 


217 


229 


2.111 


10.321 


10.241 


0.00140 


4.894 


5.391 


0.00121 


8 


366 


382 


2.087 


15.113 


14.546 


0.00000 


4.691 


5.249 


0.00405 


9 


218 


232 


2.128 


10.850 


10.915 


0.00000 


5.454 


5.856 


0.00539 


10 


201 


204 


2.030 


15.742 


15.257 


0.00166 


4.898 


5.503 


0.00491 


11 


202 


213 


2.109 


13.504 


12.891 


0.00140 


4.801 


5.217 


0.08750 


12 


25 


24 


1.920 


5.781 


5.500 


0.00000 


4.924 


5.084 


0.00000 


13 


464 


499 


2.151 


13.144 


12.703 


0.00036 


4.718 


5.390 


0.00209 



Table 2: Medium Voltage samples from the northern Netherlands Power Grid 
compared with Random graphs of the same size. 



4 Unweighted Power Grid study 

The typical study of the Power Grid as a complex system considers High Volt- 
age samples for identifying how fragile the infrastructure is. We use similar 
techniques for the Medium and Low voltage. Let us begin by recalling the basic 
complex network quantities. 

Definition 3 (Adjacency, neighbourhood and degree). If ex,y ^ E is an edge in 
graph G, then x and y are adjacent, or neighbouring, vertices, and the vertices 
X and y are incident with the edge Cx.y The set of vertices adjacent to a vertex 
X GV, called the neighbourhood ofx, is denoted by T{x). The number d{x) = 
|r(a;)| is the degree ofx. 

A global measure for a graph is given by its average distance among any two 
nodes. 

Definition 4 (Average path length (APL)). Let Vi eV be a vertex in graph G, 
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the average path length for G, Lav is: 

where d(vi,Vj) is the finite distance between Vi and Vj and N is the order of G. 

Definition 5 (Characteristic path length (CPL)). Let Vi ^ V be a vertex in 
graph G, the characteristic path length for G, Lcp is defined as the median of 
dy. where: 

is the mean of the distances connecting Vi to any other vertex Vj in G and N is 
the order of the G. 

A measure of the average 'density' of the graph is given by the clustering 
coefficient, characterizing the extent to which vertices adjacent to any vertex v 
are adjacent to each other. 

Definition 6 (Clustering coefficient (CC)). The clustering coefficient o/r„ 
zs 

^ l^(r^)l 
T„ = — -1— — 




where \E(Tv)\ is the number of edges in the neighbourhood of v and ('^2") is the 
total number of possible edges in T-u- 

This local property of a node can be extended to an entire graph by averaging 
over all nodes of the graph. 

4.1 Basic analysis 

We now consider these classic measures on the data of the Dutch Power Grid. 
We divide our data set in samples of topologically connected regions. In Table[Tl 
we report the basic analysis on the data modelled as an unweighted graphs and 
we compare each sample belonging to Low Voltage network with a random graph 
of the same size and order. The analysis for the Medium Voltage is reported in 
Table [H Referring to the table, the first column is the ID of the sample, the 
second and third represent the number of vertices N (order) and edges M (size) , 
respectively. The average degree (fourth column) is defined as <k>= The 
fifth and sixth columns report the average and characteristic path lengths, that 
is the average of the minimum distance between any two given nodes and the 
median of the same quantities, respectively. The seventh column provides an 
indication of the clustering coefficient of the nodes, that is, broadly speaking, 
an average value of the power of a node to participate in connected aggregation 
with other nodes close to it. 

We remark that the average node degree does not have highly different values 
in the Low and Medium Voltage samples, they are both around 2. Computing 
the mean over all samples' average node degree gives a value of <A:> = 2.074 
with a very small variance cr<fe> = 0.017. This value appears to be almost 
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constant considering the Low Voltage and Medium Voltage samples since the 
variance of the two categories is even smaller {(y<k>Lv = 0.016, <J<k>Mv = 
0.012). An almost constant average degree is also characteristic of the High 
Voltage Power Grid [TB|, though with a slightly higher value <k>= 2.8. This 
limited number of edges a node can manage can be regarded as a physical limit 
that each Power Grid substation has to satisfy. 

Considering path measures: Average Path Length and Characteristic Path 
Length of the Low Voltage segment of the network have generally a smaller 
path length compared to the Medium Voltage one. The clustering coefficient 
is very small especially for the Low Voltage network for which many samples 
have a zero value (i.e., absence of triangles in the graph). The difference in 
path length between the Low Voltage and Medium Voltage is due to the higher 
number of nodes the Medium Voltage network samples have while holding the 
same average node degree as the Low Voltage, together with the absence of long 
distance edges. This implies a longer path to connect any two nodes in a bigger 
network. In addition, these values of APL and GPL are in general quite high, 
if compared to other networks such as the World Wide Web. 

The clustering coefficients for the Low Voltage segment of the network are 
generally small; this is due to the strong hierarchical design of this layer of the 
physical network which resemble a tree-like structure. Gontrarily, the Medium 
Voltage segment generally presents higher values for the clustering coefficient, 
this can be justified by the different purpose the Medium Voltage network has 
in which meshed components and connection redundancies are much more likely 
to be present for robustness reasons. 

To gain an even better understanding of the tables just presented it is useful 
to compare the numbers obtained with those of Random graphs |23j and to 
identify the possible presence of Small-world properties. Small-world networks 
(SW), proposed by Watts and Strogatz in [3T], own two important aspects at 
the same time: characteristic path length close in value to the one of a ran- 
dom graph (RG) {CPLsw ~ CPLrg) but a much higher clustering coefficient 
{CCsw ^ CCrg)- Small-worlds are a better model than random graphs for 
social networks and other phenomena and thus a candidate for modelling the 
Power Grid too. To make the comparison genuine, random graphs are generated 
with the same number of nodes and edges as the real samples, imposing the re- 
sulting graphs not to have disconnected components. The values are presented 
on columns eight to ten of Tables [1] and [5] We note how the GPL of the Grid 
samples is on average twice as big as the random generated samples, thus com- 
parable to as the definition of Small- world graph according to [21] . In addition 
the clustering coefficient of the Grid samples is almost always smaller than the 
result obtained for the random generated samples; this completely contradicts 
the definition of Small- world graph according to [H] . Watts and Strogatz [3T] 
impose the following condition to the graphs they study: N ^ k ^ ln{N) ^ 1 
where N is the number of nodes, k is the number of edges per node. Such a 
condition is not satisfied by the Northern Netherlands samples and generally it 
is not satisfied by Power Grid networks as pointed by Wang et al. in [5]. In- 
terestingly, the same condition is also not satisfied by the Western States High 
Voltage Power Grid Watts and Strogatz use in [3T] and Watts analyses in [22 , 
while the results for GG and GPL satisfy the conditions for a Small-world net- 
works. Another study (i.e., [T^]) considering the European High Voltage Power 
Grid shows that the Small-world phenomenon is not shown by all the consid- 
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ered Grids, since especially the smaller (in terms of order and size) Grids fail to 
satisfy CC condition. 

In summary, the Northern Netherlands Medium Voltage and Low Voltage 
samples show a very small value of average node degree. This is mainly in- 
dependent from the size and the different purpose of the network being almost 
constant despite the different samples considered. In addition, the path length is 
quite high, given the order of the graphs, compared with other types of complex 
networks e.g., the World Wide Web. This relative high path length together 
with very small clustering properties suggests that the networks analysed do not 
strictly follow the definition of Small- world or, in terms of decentralized energy 
negotiation, it suggests that perhaps a structural change to decrease path length 
(especially the weighted one) might be necessary to empower delocalization. We 
provide an initial proposal in Section [7] on how to achieve this. 

4.2 Node Degree Distribution 

To have a general understanding of the overall characteristics of a network it is 
useful to compute certain statistical measures, one of which is the node degree 
probability distribution. More formally. 

Definition 7 (Node degree distribution). Consider the degree k of a node in a 
graph as a random variable, the function 

Nk = {veG: d{v) = k} 

is called node degree distribution. 

The shape of the distribution is a salient characteristic of the network. For 
the Power Grid, the shape is typically either exponential or a power-law. More 
precisely an exponential node degree (fc) distribution has a fast decay in the 
probability of having nodes with relative high node degree and follows a relation: 

P{k) = ae^'' 

where a and /3 are parameters of the specific network considered. While a power- 
law distribution has a slower decay with higher probability of having nodes with 
high node degree: 

P(fc) = ak"^ 

where a and 7 are parameters of the specific network considered. 

Power-law distributions are very common in many real life networks both 
created by natural processes (e.g., food- webs, protein interactions) and by artifi- 
cial ones (e.g., airline travel routes, Internet routing, telephone call graphs), pS] . 
Having a power-law distribution for node degree means that few nodes have 
a very high degree and the majority of nodes have very small degree. The 
types of networks that follow this property are referred as Scale-free networks 
([3S1[171[5S]); typical examples of Scale-free networks are the World Wide Web, 
the Internet, metabolic networks, airline routes and many others. From the dy- 
namic point of view, these networks are modelled by a preferential attachment 
model, that is, when new nodes and arcs are added to a graph, these are more 
likely to connect to nodes which have already a high degree, [551 130] ■ In addi- 
tion this network structure provides special reliability properties: high degree 
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of tolerance to random failures and, at the same time, very sensitive to targeted 
attacks towards hubs [551 1311 HH] ■ 

We compute the node degree distribution for every sample both for Low 
Voltage and Medium Voltage segments. For the most significant samples i.e., 
those belonging to the Medium Voltage and the big ones belonging to the Low 
Voltage part the node degree cumulative distribution seems to follow a power- 
law: Pk ^ . A prompt method to investigate if the node degree follows 
a power-law is to plot the cumulative node degree distribution on a log-log 
scale [32] . If the distribution in a log- log plot follows a straight line the distri- 
bution can be considered a power-law, while if the decay is faster this might 
indicate an exponential distribution. It is also possible to apply data fitting 
techniques (e.g., non- linear least square method) to identify the 7 parameter of 
a power-law. 

The most significant samples for this kind of analysis are the biggest samples 
belonging to the Medium Voltage network and the most numerous ones from the 
Low Voltage (i.e., samples #5, #9 and #10 from Tabled]). All these samples 
tend to follow a straight line in the log- log plot. Figures [2] and |3] show the dis- 
tributions for Low Voltage samples, while Figures 2] and [5] show the probability 
distribution for the Medium Voltage ones. 

We thus conclude that the Medium Voltage and Low Voltage tend to be 
Scale-free networks, although some exponential tail appear due to physical and 
economic constraints in the network. This means robustness in terms of re- 
dundancy of paths, but fragility to attacks on the hubs. The hubs tend to be 
the few nodes that most likely lead to the High Voltage segment in a certain 
geographical location. 

4.3 Betweenness 

Betweenness describes the importance of a node with respect to minimal paths 
in the graph. This is very important to identify critical components of the Power 
Grid [551 [211 m] . For a given node, betweenness, sometimes also referred as 
load, is defined as the number of shortest paths that traverse that node. 

Definition 8 (Betweenness). The betweenness 0/ vertex v & V is the number 
of shortest paths between any two vertices in graph G that contain v, i.e., 



where (Tsj(w) is the number of shortest paths from node s to node t traversing v. 

The betweenness is an important measure because it allows to find if there 
are nodes that are critical for the whole infrastructure. In fact, the removal of 
nodes with the highest betweenness can lead to critical effects on the network 
connectivity [3Tj . 

For the most significant samples in the Low Voltage network (i.e., samples 
#5 and #10) betweenness probability distribution follows an exponential decay, 
that is, the nodes with very high values of betweenness arc less likely to be 
present in the network, as shown in Figures jS] and [S] This aspect is not surprising 
since the Low Voltage network is quite hierarchical and the paths tend to follow 
the few ones admissible by the relative simple topology. In fact, the betweenness 




V 



10 



2 3 4 5 6 7 8 

K 



Figure 2: Node Degree Cumulative Probability Distribution for Low Voltage 
sample #5. Circles represent sample data, while straight line represents a power- 
law with 7 = 2.402. 




Figure 3: Node Degree Cumulative Probability Distribution for Low Voltage 
sample #10. Circles represent sample data, while straight line represents a 
power-law with 7 = 1.494. 
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Figure 4: Node Degree Cumulative Probability Distribution for Medium Voltage 
sample #2. Circles represent sample data, while straight line represents a power- 
law with 7 = 1.977. 




Figure 5: Node Degree Cumulative Probability Distribution for Medium Voltage 
sample #3. Circles represent sample data, while straight line represents a power- 
law with 7 = 2.282. 
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Figure 6: Betweenness Cumulative Probability Distribution for Low Voltage 
sample #5 (logarithmic scale). Circles represent sample data, while continuous 
line represents an exponential decay y = 0.7699e~^-^^^'^° ^ . 
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Figure 7: Betweenness Cumulative Probability Distribution for Low Volt- 
age sample #5 (semi-logarithmic scale). Circles represent sample data, while 
straight line represents an exponential decay y ~ 0.7699e~^-^^^'^° ^ . 
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Figure 8: Bctweenness Cumulative Probability Distribution for Low Voltage 
sample #10 (logarithmic scale). Circles represent sample data, while continuous 
line represents an exponential decay y = 0.6825e~^'^^®'^° ^ . 
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Figure 9: Bctweenness Cumulative Probability Distribution for Low Voltage 
sample #10 (semi-logarithmic scale). Circles represent sample data, while 
straight line represents an exponential decay y = 0.6825e~^-^^*'^° ^ . 
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Figure 10: Betweenness Cumulative Probability Distribution for Medium Volt- 
age sample #2 (logarithmic scale). Circles represent sample data, while straight 
line represents a power-law with 7 = 1.178. 
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Figure 11: Betweenness Cumulative Probability Distribution for Medium Volt- 
age sample #3 (logarithmic scale). Circles represent sample data, while straight 
line represents a power-law with 7 = 1.075. 
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probability distribution charts (Figures [S] and [5]) do not fit a straight Hue in a 
logarithmic plot, exhibiting a fast decay. This impression is also reinforced 
by the charts in Figures [7] and [HI that is, the same betweenness probability 
distribution, but on a log-linear scale: the straight line in this kind of diagram 
is a sign of an exponential decay, [32|. In addition a fitting procedure, using 
the non-linear least square method gives very good results approximating the 
betweenness probability distribution samples with an exponential function. On 
the other hand, betweenness in Medium Voltage segment seems to follow a 
power-law decay, this is shown in the logarithmic chart in Figures [TU] and 1111 
The samples from Medium Voltage network show a distribution of betweenness 
with a much fatter tail than the Low Voltage ones, that is there are nodes 
that are central in many paths. This is due to the more meshed structure the 
Medium Voltage network has, compared to the Low Voltage one. This result 
for Medium Voltage betweenness is closer to the results obtained for this same 
metric in High Voltage studies, [2H1II2- In summary, a few nodes are extremely 
critical to enable the electricity distribution to the whole network. 

4.4 Node importance 

To have a general understanding of the critical elements of the Grid, we re- 
sort to the matrix representation of the graph and study the eigenvalues and 
eigenvectors to find the most important and critical nodes of the network. 

Definition 9 (Adjacency matrix). The adjacency matrix A ~ A{G) = (ai^j) of 
a graph G is the n x n matrix given by 



Centrality refers to the importance of the node in terms of degree, between- 
ness, closeness or eigenvectors [33l|34]. In this work, wc use eigenvector central- 
ity to stress the dependence of the centrality of one node with the centrality of 
the other nodes it is connected to. The components of the dominant eigenvector 
then represent the nodes of the graph. The highest the value, the highest the 
centrality of that node, the highest the importance of that node in the Grid. 

The most critical nodes for the most interesting samples from Low Voltage 
are shown in Tables [3] and |4] while the results for samples from Medium Voltage 
are shown in Tables [5] and [B] The first column of each table represents the 
ranking position for the first ten node whose identifier is given in the second 
column. As mentioned above there are several measure to identify the most 
important nodes in a network. The core aspect of eigenvector centrality is that 
the importance of a node is dependent from the importance of their neighbouring 
nodes, and therefore to some extent from all nodes in a connected network. This 
is not the case for other measures of centrality (e.g., node degree centrality and 
betweenness centrality) whose values are not influenced by the properties of 
other nodes. 

4.5 Fault tolerance 

A related study is to evaluate the reliability of a network by analysing its con- 
nectivity when nodes are removed. There are basically two ways to perform this 




otherwise. 
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Eigenvector Centrality 


Node ID 


Ranking # 




1 


10 


2 


93 


3 


111 


4 


148 


5 


transformer 5 


6 


22 


7 


28 


8 


27 


9 


26 


10 


25 



Table 3: Eigenvector centrality ranking for Low Voltage sample #5. 



Eigenvector Centrality 


Node ID 


Ranking # 




1 


3 


2 


44 


3 


108 


4 


39 


5 


109 


6 


110 


7 


102 


8 


107 


9 


2 


10 


61 



Table 4: Eigenvector centrality ranking for Low Voltage sample #10. 



Eigenvector Centrality 


Node ID 


Ranking # 




1 


546 


2 


574 


3 


608 


4 


609 


5 


582 


6 


32 


7 


580 


8 


56 


9 


9 


10 


765 



Table 5: Eigenvector centrality ranking for Medium Voltage sample #2. 



17 



analysis: choosing the nodes randomly or selecting the nodes to be removed fol- 
lowing a certain property or metric significant for the network. Similar studies 
concerning the resilience of the High Voltage Power Grid exist, e.g., ^JM 
We apply such technique to the Medium Voltage and Low Voltage ends of the 
Power Grid using three policies for node removal: random, highest degree and 
highest betweenness driven choices. The measure that is taken into account is 
the order of the largest connected component of the network (i.e., the number of 
nodes composing the biggest cluster in the network) computed as a fraction of 
the original order of the network, and its evolution while nodes of the network 
arc removed, again the latter are considered as a fraction of the original order 
of the network. 

The random remoiia/ simulates casual errors. As shown in |35j . networks that 
follow a power-law whose characteristic parameter 7 < 3 tend to have a high 
value for the transition threshold at which they disrupt. In the samples analysed 
it seems that this is true especially for the small samples that generally have a 
cluster that is 10% of the original when almost 90% of the nodes are removed, 
as shown in Figure [T^l The situation is different for samples with higher order 
that show a cluster that is reduced to 10% of the original when about 40% of 
the nodes are removed, as shown in Figure [T^l Even if the degree distributions 
found for samples following a power-law have a parameter 7 < 3 the samples 
show a threshold effect that is more similar, according to [35], to networks whose 
characteristic 7 > 3. 

The situation is radically different when "targeted attacks" are considered. 
In particular two kind of attack policies are investigated: node degree-based 
removal and betweenness-based removal. The main difference compared to the 
random-based removal is the presence of very sharp falls that appear when 
certain nodes are targeted. The removal of selected nodes can cause a drop 
in the size of the maximal connected component even of 40%, as shown in 
Figure 1141 Node degree-based removal is much more critical than the random 
removal: by just removing 10% of the most connected nodes one reduces the 
network to only 10% of its original size. The same applies for the biggest samples 
considered both in the Low Voltage and Medium Voltage network, as shown in 
Figures [T7] and HH] for the Low Voltage and Figure [15] and [TO] for the Medium 
voltage. 

The removal of nodes based on the highest betweenness shows generally the 



Eigenvector Centrality 


Node ID 


Ranking # 




1 


351 


2 


263 


3 


324 


4 


6 


5 


12 


6 


350 


7 


299 


8 


11 


9 


80 


10 


355 



Table 6: Eigenvector centrality ranking for Medium Voltage sample #3. 
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Figure 12: Resilience for random-based removal for Low Voltage sample #2. 
The horizontal axis represents the fraction / of the nodes removed from the 
original sample; the vertical axis represents the size of the largest connected 
component 5" relative to the initial size of the graph. 
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Figure 13: Resilience for node random-based removal for Medium Voltage sam- 
ple #2. The horizontal axis represents the fraction / of the nodes removed from 
the original sample; the vertical axis represents the size of the largest connected 
component S relative to the initial size of the graph. 
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Figure 14: Resilience for node degree-based removal for Medium Voltage sample 
#10. The horizontal axis represents the fraction / of the nodes removed from 
the original sample; the vertical axis represents the size of the largest connected 
component 5" relative to the initial size of the graph. 




Figure 15: Resilience for node degree-based removal for Medium Voltage sample 
#2. The horizontal axis represents the fraction / of the nodes removed from 
the original sample; the vertical axis represents the size of the largest connected 
component S relative to the initial size of the graph. 
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Figure 16: Resilience for node degree-based removal for Medium Voltage sample 
#3. The horizontal axis represents the fraction / of the nodes removed from 
the original sample; the vertical axis represents the size of the largest connected 
component S relative to the initial size of the graph. 




Figure 17: Resilience for node degree-based removal for Low Voltage sample 
#5. The horizontal axis represents the fraction / of the nodes removed from 
the original sample; the vertical axis represents the size of the largest connected 
component S relative to the initial size of the graph. 
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Figure 18: Resilience for node degree-based removal for Low Voltage sample 
#10. The horizontal axis represents the fraction / of the nodes removed from 
the original sample; the vertical axis represents the size of the largest connected 
component S relative to the initial size of the graph. 



same behaviour, as degree-based removal, with network disruption that appear 
much faster than random-based network failures. Considering the general cor- 
relation between nodes with a certain degree and their betweenness it is not 
surprising that the two removal policies have very similar results and shape. 
The only remark that generally differentiates the betweenness-based removal 
is a little higher order of the maximal connected component compared to the 
one obtained with a degree-based removal when the same fraction of nodes is 
removed. In addition the decrease of the order of the maximal connected compo- 
nent tends to be slightly smoother than the degree-based one. Figures [TOlandl^ 
show the comparison of the two removal policies for the samples that show some 
interesting deviations in the correlation of the degree and betweenness. 

Another way to investigate the resilience is to determine which and how many 
edges have to be removed to break the graph into two disconnected components 
almost of the same order. This information is essential to have a general clue of 
the robustness of the Grid, that is how many transmission lines must be broken 
simultaneously to split the network in two components with almost the same 
amount of nodes, [19]. There are several methods and algorithms to solve this 
kind of problem that can be interpreted as a classic max-flow/min-cut problem, 
however an interesting method exploits the spectrum of the Laplacian matrix 
of a graph. 

Definition 10 (Laplacian matrix). Let D = (Dy) be a diagonal matrix with 
Da ~ d{vi) the degree of vertex Vi in graph G and A the adjacency matrix of G. 
The matrix L ^ D — A is the Laplacian matrix of graph G. 

In particular this method exploits the second smallest eigenvalue of the 
Laplacian matrix of the corresponding graph and computes the correspond- 
ing eigenvector (for a complete explanation of eigenvalues of Laplacian matrix 
spectral properties we refer to [36j[37]). This resulting eigenvector has com- 
ponents, each one representing a node of the graph, that can be cither positive 
or negative. Each node whose eigenvector component is positive belongs to one 
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Figure 19: Resilience for node degree-based and betweenness-based removal 
for Medium Voltage sample #10. The horizontal axis represents the fraction 
/ of the nodes removed from the original sample; the vertical axis represents 
the size of the largest connected component S relative to the initial size of 
the graph. Red diamonds represent the betweenness-based removal, while blue 
circles represent the node node degree-based. 




Figure 20: Resilience for node degree-based and betweenness-based removal 
for Medium Voltage sample #13. The horizontal axis represents the fraction 
/ of the nodes removed from the original sample; the vertical axis represents 
the size of the largest connected component S relative to the initial size of 
the graph. Red diamonds represent the betweenness-based removal, while blue 
circles represent the node node degree-based. 



23 



sub-graph, while the ones with negative components belong to the other sub- 
graph. The edges that connect nodes belonging to the different sub-graphs are 
the critical edges that if removed lead to two different sub-graphs. The more 
edges connect the two sub-graphs, the more robust is the Grid. Once the two 
sub-graphs are identified it is possible to iterate the method on each sub-graph 
and find again the most critical edges. 

The number of critical edges that, if removed at the same time, disrupt the 
network evenly in two or more sub- networks, are shown in Table [7] and Table [5] 
for Low Voltage and Medium Voltage respectively. The first column represent 
the sample identifier, while the second column represent the number of edges 
to be removed simultaneously to split evenly the network. Generally Medium 
Voltage networks are more robust to edge failures than Low Voltage with the 
number of edges to be removed that is double. Medium Voltage sample #2 is 
the most reliable to line disruption attack and indeed this is an indication of 
the high mashed structure the biggest sample analysed owns. 



Sample # 


Number of Critical 




Edges 


1 


2 


2 


2 


3 


2 


4 


1 


5 


2 


6 


2 


7 


1 


8 


1 


9 


2 


10 


3 


11 


1 



Table 7: Number of critical edges according to Laplacian eigenvalue method for 
Low Voltage samples. 



Sample # 


Number of Critical 




Edges 


1 


2 


2 


27 


3 


4 


4 


5 


5 


3 


6 


4 


7 


4 


8 


1 


9 


6 


10 


4 


11 


4 


12 


1 


13 


6 



Table 8: Number of critical edges according to Laplacian eigenvalue method for 
Medium Voltage samples. 

In summary, the results for the Low Voltage and Medium Voltage show 
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disruption behaviours. These networks are quite immune to random failures to 
which the networks present a constant degrading disruption, while they deeply 
suffer from certain characteristic nodes to be removed. 

5 Weighted Power Grid study 

The purely topological study of the Power Grid just presented already gives im- 
portant information about the connectivity and robustness of the Medium and 
Low Voltage Grids, though it docs not consider the different physical properties 
of the cables. These can vary greatly for different sections of the Grid and pro- 
vide essential indications to establish the behaviour of a link. Next we perform 
an analysis of the same samples of the Grid also considering the resistance as 
the weight of the graph model of Definition [2j 



ID 


Weighted 
Characteristic 
Path Length 


Edge Average 
Weight 


Normalized 
Weighted 
Characteristic 
Path Length 


1 


2.000 


0.698 


2.865 


2 


1.429 


0.595 


2.402 


3 


3.066 


0.739 


4.149 


4 


3.087 


0.699 


4.414 


5 


12.136 


0.741 


16.378 


6 


3.889 


1.648 


2.360 


7 


4.162 


0.348 


11.960 


8 


5.112 


0.876 


5.836 


9 


7.872 


0.583 


13.503 


10 


6.407 


0.785 


8.162 


11 


2.967 


0.592 


5.012 



Table 9: Weighted analysis of the Low Voltage samples from the northern 
Netherlands Power Grid. 

Take the samples analyzed in Tables [T] and [2l but now consider the weighted 
graph definition. The notion of a characteristic path length can be extended to 
take the weights into account yielding the values shown in Tables ^ and [TOl In 
each table, the second column contains the characteristic path length resulting 
in the weighted graph (WGPL), formally: 

Definition 11 (Weighted characteristic path length (WGPL)). The weighted 
characteristic path length for graph G, L.^cpi is the median for all {vi,Vj) G V 
of the following distance 

dw{pi: ) — ^ ^ ^Ws^t 

such that Cw, t is an edge in the minimal weighted path between Vi and Vj. 

The third column provides the average value of the weights of all edges; while 
the fourth column shows a normalized value for the weighted characteristic path 
length (NWGPL) obtained by dividing the WGPL by the average weight of the 
edge belonging to the same data sample. This normalization is performed to 
have a measure to compare the unweighted and the weighted samples whose 
results are shown in Section [51 
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ID 


Weighted 
Characteristic 
Path Length 


Edge Average 
Weight 


NormaUzed 
Weighted 
Characteristic 
j^axn ijerigTin 


1 


185.916 


12.779 


14.549 


2 


108.011 


11.851 


9.987 


3 


153.402 


8.608 


17.821 


4 


163.067 


9.217 


17.692 


5 


127.258 


7.122 


17.868 


6 


134.661 


13.106 


10.275 


7 


187.084 


16.382 


11.420 


8 


148.058 


7.193 


20.584 


9 


99.385 


7.421 


13.392 


10 


126.845 


6.850 


18.518 


11 


92.060 


8.764 


10.504 


12 


38.084 


6.915 


5.507 


13 


232.475 


13.810 


16.834 



Table 10: Weighted analysis of the Medium Voltage samples from the northern 
Netherlands Power Grid. 

Due to the relative short length of the Low Voltage networks cables, the 
WCPLs for this segment of the network are small, as well as the average weight 
of each edge (almost all of them are below the unit). The situation is different 
for the Medium Voltage networks which are higher since the cables and paths 
span across wider geographical areas. The discrepancy can be explained by 
the different purpose for which these networks are designed: a bridge network 
from High Voltage transmission lines and end-user distribution (Medium Volt- 
age network) and the final end delivery (Low Voltage network). In fact, both 
the WCPL and the edge average weight for Medium Voltage samples are ap- 
proximately two order of magnitude greater than the Low Voltage ones. This is 
indeed due to an extension of Medium Voltage cables that range from hundred 
meters to kilometres, while Low Voltage extend usually around tens of metres. 

5.1 Weighted Node Degree Distribution 

Though no value is associated to a node, the weights of the incident edges also 
influence the node properties. One way of seeing this, is by defining a weighted 
node degree. 

Definition 12 (Weighted degree). Let x £ V be a vertex in a weighted graph 
G, the weighted degree of x, dw{x) is: 

dw{x) = ^ 'Wx,y 

yeV(x) 

where Wx,y is the weight of the edge joining vertices x and y and T{x) is the 
neighbourhood of x. 

The weighted distribution is straightforwardly obtained by using the weighted 
degree in Definition [T] 

For the most significant sample of the Low Voltage, as shown in Figures 
and l221 the shape of the distribution is close to an exponential one with a quite 
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Figure 21: Weighted Node Degree Cumulative Probability Distribution for Low 

Voltage sample #5 (logarithmic scale). Circles represent sample data, while 

continuous line represents a sum of exponential decays y = O.SOTSe""'^^^^^ + 
0.0904e-°-i37to_ 

fast decay. The situation looks different in Medium Voltage samples. The very 
first part of the distribution is well fitted by an exponential shape, while the 
central part of the distribution, and especially the tail, fit best a power-law like 
shape as visible in Figures [23] and An explanation of such behaviour between 
the the most numerous samples of the two ends of the Grid is due to the order 
and size of the Medium Voltage samples which are from two to four times bigger 
than the Low Voltage samples, thus having a higher likelihood of far different 
values in weighted node degree. 

5.2 Betweenness 

Considering the weighted definition of path it is possible to compute between- 
ness in the weighted scenario and again betweenness of a node can be seen as 
a random variable thus obtaining the corresponding probability distribution. 
The shape of the distribution does not change much compared to the same un- 
weighted samples for the Low Voltage network, as shown for samples #5 and 
^10 in Figures 1251 and [26l the distribution is best approximated by an expo- 
nential decay or by a sum of exponential contributions. For Medium Voltage 
samples, the changes between unweighted and weighted paths influence the be- 
tweenness probability distribution whose shape in these conditions seems to be 
better approximated by an exponential or sum of exponential components as 
shown in Figure 1271 This change in the distribution of the number of shortest 
paths that traverse a node between the weighted and the unweighted graph is 
clearly an indication that some property change between the two analysis and 
it is worth to remember that the weighted path analysis better approximate the 
actual routes current flows follow. 



27 




Figure 22: Weighted Node Degree Cumulative Probability Distribution for Low 
Voltage sample #10 (logarithmic scale). Circles represent sample data, while 
continuous line represents a sum of exponential decays y = 0.1538e~^^-^'''^ + 
0.8378e-o-4909^ 




Figure 23: Weighted Node Degree Cumulative Probability Distribution for 
Medium Voltage sample #2 (logarithmic scale). Circles represent sample data, 
while straight line represents a power-law with 7 = 1.354. 
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Figure 24: Weighted Node Degree Cumulative Probability Distribution for 
Medium Voltage sample #3 (logarithmic scale). Circles represent sample data, 
while straight line represents a power-law with 7 = 1.374. 




Figure 25: Betweenness Cumulative Probability Distribution for Low Voltage 
sample #5 considering the weighted graph (logarithmic scale) . Circles represent 
sample data, while continuous line represents a sum of exponential decays y = 
-0.1051e3-38iio^°^ + 0.8084e-i-3i7io-*a;^ 
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Figure 26: Betweenness Cumulative Probability Distribution for Low Volt- 
age sample #10 considering the weighted graph (logarithmic scale). Circles 
represent sample data, while continuous line represents an exponential decay 
y = 0.6456e-6-"'3-i°"'^. 
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5.3 Node Importance 



The criticality or importance of nodes in the network is best studied using 
weights. In fact, an edge with high capacity (i.e., weight) makes a node very 
important, conversely many edges with httle capacity (i.e., weight) make a node 
almost irrelevant. The approach considered is similar to the one leading to 
the eigenvector centrality computation performed in Section 14.41 However, to 
compute this metric, a weighted form of the adjacency matrix is necessary (cf. 
Newman [55]). 

Definition 13 (Weighted adjacency matrix). The weighted adjacency matrix 
Am ~ Au]{G) ~ of a graph G is the n x n matrix given by 

{Wij if ViVj € E and has weight Wij, 
otherwise. 

Computing the eigenvector corresponding to the principal eigenvalue, one 
obtains a ranking among the various nodes of the network. The results for 
eigenvector centrality for two Medium Voltage samples are shown in Tables [11] 
and [T^] while for two Medium Voltage samples arc shown in Tables [T2] and [TJ] 
The first column of each table represents the ranking position for the first ten 
node whose identifier is given in the second column. It is interesting to note 
that with this ranking the first ten most important nodes of the graph, and to 
a certain extent critical substations for the Power Grid, are in general differ- 
ent than those of the unweighted analysis. This aspect reinforces a consistent 
difference in node properties between the two type of analysis performed. 



Eigenvector Centrality 
Ranking # 


Node ID 


1 


192 


2 


191 


3 


6 


4 


24 


5 


137 


6 


130 


7 


135 


8 


129 


9 


178 


10 


transformer 3 



Table 11: Eigenvector centrality ranking for Low Voltage sample ^5 considering 
the corresponding weighted graph. 



5.4 Fault tolerance 

Fault tolerance can be evaluated based on the removal of nodes following strate- 
gies similar to the unweighted case. Since the random removal yields exactly 
the same result for the weighted and unweigted case, here we focus on the node 
degree-based removal policy which considers the weighted node degree defini- 
tion. The disruption behaviour of the network samples is very similar to the 
unweighted situation: the network suffers deeply these targeted attacks; a very 
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Eigenvector Centrality 


Node ID 


Ranking # 




1 


3 


2 


44 


3 


108 


4 


5 


5 


107 


6 


55 


7 


110 


8 


109 


9 


102 


10 


39 



Table 12: Eigenvector centrality ranking for Low Voltage sample #10 consider- 
ing the corresponding weighted graph. 



Eigenvector Centrality 
Ranking # 


Node ID 


1 


577 


2 


580 


3 


575 


4 


822 


5 


702 


6 


706 


7 


546 


8 


578 


9 


170 


10 


90 



Table 13: Eigenvector centrality ranking for Medium Voltage sample #2 con- 
sidering the corresponding weighted graph. 



Eigenvector Centrality 


Node ID 


Ranking # 




1 


324 


2 


351 


3 


263 


4 


350 


5 


299 


6 


410 


7 


393 


8 


48 


9 


124 


10 


80 



Table 14: Eigenvector centrality ranking for Medium Voltage sample #3 con- 
sidering the corresponding weighted graph. 
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Figure 28: Resilience for weighted node degree-based removal for Low Voltage 
sample ^^5. The horizontal axis represents the fraction / of the nodes removed 
from the original sample; the vertical axis represents the size of the largest 
connected component S relative to the initial size of the graph. 



small percentage of removed nodes causes an important loss in the size of the 
biggest component left in the network. The comparison between Figures [551 
and [17] for the Low Voltage samples, and Figures [29] and [15] for the Medium 
Voltage samples provides a general correlation between high degree nodes in 
the unweighted graph and high degree nodes in the weighted one. If one takes a 
closer look at the disruption charts for the same samples some small differences 
can anyway be noticed. The nodes with the highest weighted degree cause a 
bigger damage to the network when removed in the very first iteration than 
nodes with higher degree in unweighted networks, this behaviour is shown in 
Figures [30] and [21] The situation then changes in the later stages of the removal 
process when a bigger disruption is caused by nodes with higher node degree in 
traditional sense. 



6 Unweighted vs. Weighted comparison study 

The weighted study of the Power Grid presented in the previous section has 
already highlighted the more precise information available with such study. Next 
we consider more in detail the comparison between the unweighted and weighted 
study for the most indicative measures. 

Considering a minimal path, one may wonder if introducing weights changes 
the number of traversed nodes by such minimal paths on average. Even if 
the WCPL and the NWCPL values are close, there might be differences in 
the number of nodes and (therefore edges) traversed when following a minimal 
path. This is particularly interesting from the practical point of view, as it 
indicates the number of transformers and distribution substations traversed in 
the Power Grid. These points are critical in terms of additional losses that 
are associated with substations and transformers, and in turn in the number 
of potential points of failure that a path traverses. Figures [35] and [33] show 
the results for the Low Voltage and Medium Voltage, respectively. Each bar 
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Figure 29: Resilience for weighted node degree-based removal for Medium Volt- 
age sample jj^2. The horizontal axis represents the fraction / of the nodes 
removed from the original sample; the vertical axis represents the size of the 
largest connected component S relative to the initial size of the graph. 




Figure 30: Resilience for node degree-based removal for Medium Voltage sample 
#3. The horizontal axis represents the fraction / of the nodes removed from 
the original sample; the vertical axis represents the size of the largest connected 
component S relative to the initial size of the graph. Red diamonds represent 
the weighted node degree-based removal, while blue circles represent traditional 
node degree-based removal. 
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Figure 31: Resilience for node degree-based removal for Low Voltage sample 
#10. The horizontal axis represents the fraction / of the nodes removed from 
the original sample; the vertical axis represents the size of the largest connected 
component S relative to the initial size of the graph. Red diamonds represent 
the weighted node degree-based removal, while blue circles represent traditional 
node degree-based removal. 
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Figure 32: Percentage increase in number of node traversed between weighted 
and un- weighted Low Voltage samples. 
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Figure 33: Percentage increase in number of node traversed between weighted 
and un- weighted Medium Voltage samples. 



represents the average percentage increase in the number of nodes traversed 
along the shortest path between any two nodes for the unweighted and the 
weighted situation. It is interesting to note that for several samples of the Low 
Voltage there is no difference in the number of traversed nodes, thus reinforcing 
the idea of a highly hierarchical tree-like structure whose paths are fixed by 
the built-in topology of the Grid independently of the associated edge weights. 
The situation though is quite different for the Medium Voltage. In fact there is 
an increment of traversed nodes between the weighted and unweighted models 
(especially for the meaningful samples) on average of about 50%. This is a 
clear indication of a meshed network for which there are less imposed paths 
and in which the weights have an important role. It is important to notice that 
for the biggest sample (almost 900 nodes) the number of visited nodes when 
following a path between any two nodes increases by more than 80% comparing 
the unweighted and the weighted situation. 

Considering the node degree distribution it seems that the weighted analy- 
sis tends to reduce the contribution of the tail components of the distribution, 
thus being more compact especially for Low Voltage samples. This is due to 
the small variance of the weighted node degree that especially these low-end 
samples of the Grid show. Samples that in the unweighted analysis show a 
power-law distribution, when considered weighted tend to assume an exponen- 
tial form or a sum of exponential contributions. The biggest samples analysed 
in the Medium Voltage (samples #2 and #3) tend to even have a faster decay 
than the unweighted situation. The same consideration applies to the between- 
ness probability distribution which for the Low Voltage samples still shows a 
best approximation by an exponential decay. For the samples belonging to the 
Medium Voltage the same tendency appears: a more compact distribution of 
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the number of paths traversing nodes. 

Considering a representation of the relationships between nodes and their 
weights, the ranking of importance obtained with the eigenvector centrahty 
computation between nodes changes substantially. This is understandable since 
the weight deeply influences the properties, and therefore importance, of a node. 
This dissimilarity in node characteristics is also shown by the different behavior 
the removal of the highest degree nodes brings to the network connectivity. 
Although the general behavior functioning similar, the first nodes removed with 
the highest weighted degree tend to have a more damaging impact than the 
corresponding policy removal in an unweighted graph, this behavior then tends 
to reverse while removing more and more nodes: the most damaging results in 
term of network connectivity is then caused by nodes with highest degree in the 
unweighted network. The comparison between weighted and unweighted node 
degree-based removal is shown in Figure [21] for a Low Voltage sample and in 
Figure 1201 for a Medium Voltage sample. 

7 Topological influence on Energy Exchanges 

Traditionally, energy has been 'pushed' hierarchically from the large scale pro- 
duction facilities to the end users. Famous is the quote of Samuel Insull (XIX 
century) : 

"every home, every factory and every transportation line will obtain 
its energy from one eommon souree, for the simple reason that that 
will be the cheapest way to produce and distribute it. " 

Clearly, things are rapidly changing on today's Power Grid and new models are 
emerging where delocalized production is the norm, rather than the exception. 
The trend will call for an infrastructure that supports energy trading among 
any actor connected to the Power Grid. 

The Complex Network Analysis that we provided so far gives a statistical 
aggregated view of the current infrastructure for the Low and Medium Volt- 
age Grids. The natural next question that arises concerns the usability of such 
infrastructure for the delocalized energy exchange. To answer it, we propose 
to tie statistical properties of the Power Grid with a trading cost. The cost 
represents a balance frontier below which the actors are motivated to trade and 
above which they are not. It is important to remark that we do not claim of 
having identified "the" cost function, but rather we propose that Complex Net- 
work Analysis measures are deeply connected to the success of a decentralized 
energy market. 

7.1 Relating topology to price 

In general, establishing energy pricing is not a simple task since several aspects 
influence the overall price at which electricity is sold. There are aspects con- 
nected to the supply side such as fuel prices, policy regulations, load losses and 
bidding strategies; on the other hand, there are elements connected to the de- 
mand side such as human behaviors, natural phenomena that influence habits 
and thus consumption. Recent proposals and methods for price allocation in- 
clude nodal pricing |39] . which is particularly indicated for distributed genera- 
tion solutions because of the price benefits it brings to the customer [30] . It is 
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also important to notice that the savings deriving from distribution losses can 
be extremely important [HJ |3H US] . A set of factors is most definitely tied to 
infrastructural properties of the distribution network, as illustrated for instance 
in the economic studies of Harris and Munasinghe [HI SS] , most notably: 

• losses both in line and at transformer stations, 

• security and capacity factors, 

• line redundancy, and 

• power transfer limits. 

The listed technical parameters arc naturally associated with a topological pa- 
rameter, namely: 

• Line losses are related to and thus expressed as a function of the weighted 
characteristic path length Lune^ = HWCPLn) 

• Substation losses are expressed as a function of the (average) number of 
nodes visited while travelling from source node to destination node along 
the weighted shortest path Lsubstationij = f{\WSPij\). The significant 

dimension is Lsubstationr, = JUVVSPnI) where \WSPn\ is the average 
number of nodes traversed in a shortest path in the network N. 

• Robustness from a topological perspective is expressed as the fraction of 
the maximal connected component compared to its original size once a 
certain fraction of nodes is removed. RohM = f{N,p) where N is the 
network under evaluation and p is the removal policy adopted. 

• Line redundancy is simply mapped to a topological metric that counts 
the number of paths (without cycles) that are available between any two 
nodes and the cost associated to this redundancy Redij = f{\Pij\,Wij) 
where is the set of paths between nodes i and j and Wij is the weight of 
the worst case redundant path. A global metric for network N is RccLn = 

/(p^,u^),Vi,je 

• Network capacity may be considered as a function of the weighted charac- 
teristic path length where the weight is the maximal supported operating 
current of the cable CapN ~ fiWCPLpf) 

These topological ingredients provide two sorts of measures, the first one a 
gives an average of the dissipation in the transmission between two nodes 

— f {P line N ^ substation (1) 

the second one /? is a measure of reliability/redundancy in the paths among any 
two nodes 

13 = f{RobN,RedN,CapN). (2) 

We argue that these two factors influence the inclination of prosumcrs (energy 
consumers/producers) to trade energy on the Power Grid. In fact a high value 
of the a parameter represents a high level of losses experienced for transporting 
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Transport cost of energy based on topologcal properties for Low Voltage 




Figure 34: Transport cost of energy based on the topological properties for Low 
Voltage samples. 
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Figure 35: Transport cost of energy based on the topological properties for 
Medium Voltage samples. 
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energy in the network, either in distribution hnes or substations. Additionally, 
the reliability and ability to bring sufficient energy to the end users plays an 
important role. In fact, if proper levels of robustness of the network or resilience 
to failures are not the norm, the prosumer inclination to sell energy as well as the 
end user to buy it will be limited. Furthermore, if the availability of redundant 
paths for electricity routing in case of partial disruption of the network arc 
insufficient leads to a high value of /3, and in turn a disincentive for trading. 
To better understand the constituents of a and (3, we consider next a possible 
instantiation of these parameters using the data of the Dutch Power Grid. 

7.2 An example 

To give an impression of how the parameters can be used to assess the success of 
the energy market, we provide an example next. We stress that this is simply an 
example, and it does not have the ambition to provide the successful parameter 
for the delocalization of the energy distribution market. 

• Losses on the transmission/distribution line can be expressed by the quo- 
tient of the weighted characteristic path length and the average weight of 
a line (a weighted edge in the graph): 

T _ WCPLt 



• Losses at substation level are expressed as the number of nodes (on av- 
erage) that arc traversed when computing the weighted shortest path be- 
tween all the nodes in the network; 



LsubstatioUN = NodeswcPLN 

Robustness is evaluated with random removal strategy and the weighted 
node degree-based removal by computing the average of the order of max- 
imal connected component between the two situations when the 20% of 
the nodes of the original graph are removed. It can be written as: 

Redundancy is evaluated by covering a random sample of the nodes in 
the network (40% of the nodes whose half represents source nodes and the 
other half represents destination nodes) and computing for each source 
and destination pair the first ten shortest paths of increasing length. If 
there are less than ten paths available, the worst case path between the two 
nodes is considered. To have a measure of how these resilient paths have 
an increment in transportation cost, a normalization with the weighted 
characteristic path length is performed. We formalized it as: 



Red 



^ ~ WCPL 



Network capacity is considered as the value of the weighted characteristic 
path length, whose weights are the maximal operating current supported, 
normalized by the average weight of the edges in the network (average 
current supported by a line). That is: 
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With these instantiations, equations ([T]) and ([2]) become: 



1 ^substations ) — ^l' 



'lincN ^~ -^substation 

Redjsf 



■N 



(3) 



/3 = f{RedN,RohN,Capisi) 



RobN ■ ln{CapN) 



(4) 



The functions to compute a and /3 are only few of many available possible 
ones. The choice made here is to have a simple mechanism to assess the potential 
exchange costs of different networks. Equation ([3]) is a basic sum over the losses 
that are experienced both at line and at substation level. Equation (jj]) takes 
into account the aspects of reliability and tolerance of the network: the higher 
the (3 the more prone to failures and less reliable the network is. The cost of 
increasing the number of paths to provide more redundancy is the dividend in 
the fraction, while elements improving reliability act as divisor. 

With these quantities, one can form an impression of what the influence of 
the cost of transportation is for the decentralized energy exchange. If the cost is 
too high because of an infrastructure with high chances of failure (high /3) and 
high resistance (high a), then the decentralized market will not be incentivized. 
On the other hand, for low a, /3, it will be economically attractive to have a 
decentralized energy market. In Figures l34l and 1551 we show a combination of 
a, fj obtained with the functions described above and with an hypothesis of a 
quadratic increase of energy price with the increment in a and /?. We report the 
position of the analysed samples as white circles in FigureslMlandlSSlrespectivelv 
for Low Voltage and Medium Voltage samples. By performing an economic 
study, which we stress is beyond the scope of the present treatment, one then 
can identify what the threshold is for the feasibility of a decentralized market 
(grey area in Figure [36]) and then conclude what topological modifications are 
necessary to the Medium Voltage and Low Voltage infrastructure in order to 
allow the energy exchange. 



Complex Network Analysis studies are becoming more and more popular given 
the amount of natural and human based complex systems. The Power Grid is 
clearly amenable to such study and a number of these have been performed on 
the High Voltage Grid around the globe. These have been an inspiration for the 
present work which though is novel in its analysis of the Medium Voltage and 
Low Voltage, in its use of weights, and its motivation stemming from energy 
negotiation by micro/medium prosumers. Next we survey the most notable 
related work. 

Albert et al. [46] study the reliability aspects of the United States Power 
Grid. In extreme summary this work is particularly relevant for the big sample 
it analyses representing the whole North American High Voltage Power Grid 
and focusing in illustrating the cascading effects evaluating a connectivity loss 
property the authors define. The effects affecting the Power Grid based on this 
metric show important differences between the various policies of node removal 
(random, node degree or betweenness based). Betweenness computation is re- 
markable as it is used to identify the important nodes in the network, however 
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Transport cost of energy based on topological properties for Medium Voltage 




Figure 36: Transport cost of energy based on the topological properties for 
Medium Voltage with supposed economic convenience threshold (grey thick 
line) . 

the article does not take into account any sort of weight to be associated to the 
power lines. 

Crucitti et al. [T3] analyse the Italian High Voltage Power Grid from a topo- 
logical perspective. This work is particularly relevant for the concept of effi- 
ciency that is used to understand the performance of the network. This metric 
is evaluated as a function of the tolerance of load both for edges and nodes. It 
is interesting that some sort of weights are used for this analysis: a capacity 
measure is associated with nodes while weight is associated with edges based on 
residual capacity of nodes. In this article too, some strategies of failure simula- 
tion are taken into account (random and betweenness- based removal). However, 
the size of the sample analysed is small compared to other works and the type 
of weight attributed to edges is not related to any physical quantity (e.g., lines 
resistance), but is only related to topological betweenness. 

A different model is presented by Chassin et al. [12] where the analysis fo- 
cuses on the North American Power Grid. The authors start with the hypothesis 
that the Grid can be modelled as a Scale-free network. This work is extremely 
relevant for the large size of the sample analysed (more than 300000 nodes) and 
the use of reliability measures that are typical of power engineering (e.g., loss 
of load probability) to quantify the failure characteristics of the network from a 
topological point of view. The similarity of the results obtained by considering 
reliability with authors' topological measures and other non topological studies 
for electrical Grids is interesting. However, a study of betweenness is unfortu- 
nately missing. Given the size of the sample (although it is not explicitly stated 
if Medium Voltage components are considered or not) it would have useful to 
compute the betweenness of the nodes in order to understand if and how the 
betweenness behaves in this sort of big sample. 

Holmgren [T^l analyses the Nordic Power Grid involving the High Voltage 
Grids of Sweden, Finland, Norway, and a great part of Denmark comparing 
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Table 15: Comparison table between literature studies using Complex Network 
Analysis applied to Power Grid networks. 

these with the U.S. Power Grid. A resilience analysis is performed and with 
the inclusion in this work of some fictitious scenarios of failure of the Grid 
and the possible adoptable solutions together with their resulting benefits. A 
computation of the betweenness property of the graph might have been useful 
to understand the differences between the different samples; also a weighted 
graph study might have pointed out even more interesting aspects between the 
various networks, but these are missing. 

Casals et al. [17] analyse the whole High Voltage European Power Grid 
and try to extract non-topological reliability measures investigating the topo- 
logical properties of the network. The Power Grid analysed is the High Voltage 
end composed of almost 2800 nodes that span across the European continent. 
The assumption is that node degree distribution follows an exponential decay 
for every single network composing the European network, each one having a 
characteristic parameter that is related to the robustness of the specific Grid. 
Although this study is based on relevant samples, but half of the considered 
Grids are small both in size and order (below 100 nodes). The most interesting 
aspect is the use of new indicators to assess network reliability, but, as the same 
authors explicitly state, these metrics need more test and a deeper study. As re- 
marked for other works, there is no mention about using weights to characterize 
the edges in the networks. 

Casals et al. [TB] consider the High Voltage Grids of many European countries 
analyzing them together and as separate entities. This work has an overall rele- 
vant sample although no information are given for each single Grid which might 
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have smaller significance when analysed alone. The most interesting aspects 
are the evaluation of Small-world properties for the networks composing the 
European Grid and the resilience test with both random and node degree-based 
removal strategies. There is no mention about using weights to characterize the 
edges in the network and no betweenness computation to find out critical nodes 
in term of paths covered. 

Sole et al. [17] go further in exploring the same Power Grid data analysed 
in [TB], in particular they focus on targeted attacks. The model is based on the 
assumption, also verified by empirical data, that there is no correlation between 
nodes having a certain degree to be connected to each other. This work has 
an overall relevant sample and it focuses on simulating difFerent failure events 
(random or targeted), establishing an interesting correlation between topologi- 
cal and non-topological reliability studies. The major point of improvement is 
related to the small size of half of the samples used (below 100 in order) and 
the possibility of introducing weights for edges related to some of the physical 
properties. In addition, an evaluation of the betweenness is missing in order to 
understand if other nodes appear to be critical and so to be targeted using this 
different removal metric. 

Crucitti et al. [50] analyse the High Voltage Power Grid of Italy, France and 
Spain to detect the most critical lines and propose solutions to address their 
vulnerabilities. This work has some very valuable aspects such as the comparison 
of the Grids of three different countries (i.e., Italy, France and Spain), the 
identification of the most vulnerable edges, the damage provoked by an attack 
and possible improvements based on the efficiency metric. The sample used is 
considerably small; in addition there is no use of weights to characterize the 
edges. Thus it is not possible to discover which is the weight of the most critical 
edges identified and if there is a correlation with the unweighted analysis. 

Rosato et al. [T^] analyse the same network samples studied in [5^ to inves- 
tigate the main topological properties of these Grids (i.e., Italian, French and 
Spanish High Voltage Grids). The contributions of this work include the com- 
parison of the Grid of three different countries, the identification of the most 
vulnerable edges, and the damage of an attack and achievable improvements 
based on adding edges. It also studies the node degree distribution and the 
shortest path length distribution for these samples. It is interesting to note 
how the authors clearly show the correlation between country geography and 
topological measures. The sample used is notably small; in addition there is no 
use of weights to characterize the edges. 

Watts [24] dedicates a subsection of his book to exploring the properties of 
the Western States Power Grid. The study gives motivations to the Small-world 
modeling. The analysis focuses on specific metrics such as network contraction 
parameters and the comparison between different models (i.e., relational and 
dimensional models). Therefore being Small- world the focus of the analysis, 
other typical Complex Network Analysis measuring arc not performed (e.g., 
node degree distribution, betweenness distribution). 

Wang et al. [5] study the Power Grid to understand the kind of commu- 
nication system needed to support the decentralized control required by the 
Smart Grid. The analysis is based both on real Power Grid samples and syn- 
thetic reference models belonging to the IEEE literature. This work has some 
very valuable aspects such as the investigation of a significant sample, the in- 
dividuation of a new model for the node degree probability distribution, and 
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the investigation of the physical impedance distribution of the Grid samples. 
All these factors bring to the development of a new model to characterize the 
Power Grid. An aspect that might have been analytically evaluated is the path 
length in the various samples analysed and the betweenness computation to 
characterize even those distribution analytically. The use of electrical proper- 
tics is extremely interesting, however the analysis performed is dissociated to the 
physical graph properties therefore not considering a weighted graph structure. 

There are also some brief studies related to the Power Grid that appear as 
examples in more general discussions about Complex Networks. In particular, 
Amaral et al. in [32] show a study of the Southern California Power Grid and 
the model following an exponential decay for node degree distribution. Watts 
and Strogatz in |21j show the Small-world phenomenon applied to the Western 
States Power Grid while Newman, within a more general work [?5], shows the 
exponential node degree distribution for the same Grid, while Barabasi et al. |49j 
model the Power Grid as a Scale- free network characterized by a power-law node 
degree distribution. 

As shown in Table [T5] the analysis performed in the present study has a 
sample significant in size (considering the singles samples altogether) higher than 
most other studies of the High Voltage network. The new aspects that make the 
present analysis unique are the focus on the Medium Voltage and Low Voltage 
network that no other study takes into account, and also the usage of weighted 
graph that other studies miss. In addition, the analysis is complete investigating 
betweenness and resilience analysis through the computation of various failure 
policies. Nonetheless the analysis examines the Small-world property for the 
Medium Voltage and Low Voltage network showing how the samples under 
exam do not belong to this class of networks. 

9 Conclusions 

When facing a global complex system such as the Power Grid, it is necessary 
to combine precise tools to consider the local phenomena with global statistical 
tools that provide an overall view. In the present study, we look with a weighted 
model at the Medium and Low Voltage Grids with the aim of understanding its 
potentials as infrastructure for delocalizcd energy distribution. The study takes 
the North Netherlands data as basis for its analysis and proposes and adapts a 
number of statistical topological measures for the specific goal. 

A number of novelties arc a trademark of the current proposal: the study of 
the lower layers of the Power Grid, the study of energy distribution rather than 
simply resilience, the use of a weighted topological model and, most notably, the 
proposal of tying the topological properties to values denoting the attractiveness 
for the prosumers to trade energy. But let us be more specific in summarizing 
the results of the present treatment. 

An interesting aspect is the difference experienced in the shortest path be- 
tween the weighted and unweighted graph. The difference is most notable in 
the number of nodes that need to be traversed. This is interesting since the 
weighted paths considering resistance are closer to the ones really travelled by 
the energy flow, therefore traversing much more nodes than the ideal situation 
can lead to additional losses in substations and transformers together with a 
larger number of potential point of failure. 
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Considering graph statistical properties, the node degree distributions tend 
to follow a power-law (at least for the most significant samples), that is there 
are few nodes that have many connections, while the majority has a very limited 
number. As the literature shows there is quite a common consensus about the 
type of node degree distribution High Voltage Power Grids have, i.e., exponen- 
tial, while a study that takes into account many more nodes obtains a power- law 
distribution [T2] . The explanation of the Medium and Low Voltage Grids may 
reside in the relatively small number of stations and transformers that receive 
their input from the High Voltage network and have to distribute this input to 
many more substation at lower voltages. 

The betweenness plots have a specific characteristic for the Low Voltage 
samples: unlike what presented in the literature reporting a power-law charac- 
terizing the High Voltage, the Low Voltage follows an exponential decay. The 
topology itself of the Low Voltage network, in which the paths are much more 
forced due to the greater hierarchy of the Low Voltage network, implies a be- 
tweenness distribution with a more compressed tail than a fat-tailed power-law. 
While a more meshed network as the Medium Voltage has more nodes that take 
part in different number of paths. 

Another remarkable aspect is the relatively higher tolerance that is shown 
by the Medium Voltage network when edges are removed compared to the Low 
Voltage network. This is due to the more meshed structure of the Medium 
Voltage which is therefore less prone to failures than the Low Voltage. In fact, 
the average number of links that need to be eliminated at the same time to cut 
the network in two or more components are more than the ones for Low Voltage. 

From these findings, we move towards our claim, that is, the influence of 
these global topological factors on the actual possible use of the Power Grid as 
an infrastructure for delocalized energy exchange. We make an initial proposal 
of how the Complex Network Analysis topological metrics can be clustered in 
two parameters a and /3 that need to be optimized in order to facilitate energy 
negotiation among prosumers. To optimize them, one may have to change the 
topological properties of the Grid. It might turn out that to reduce the WCPL 
needs to lay more cables between rural and urbanized areas, or to put more 
interconnections between neighboring urbanized areas. 

Clearly there is much more to be investigated in the direction of Complex 
Network Analysis for the Power Grid. The models need to be enriched with 
ever more Power Grid specific characteristics. The dynamics of the growth of 
the Power Grid may also provide new insights. Simulations of possible future 
scenarios will also help to identify the right infrastructure of the future (Smart) 
Grid, not to mention the importance of coupling such technical studies with 
their economic counterparts. 
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